Spin chains with electrons in Penning traps 

G. Ciaramicoli, I. Marzoli, and P. Tombesi 

Dipartimento di Fisica, Universitd degli Studi di Camerino, 62032 Camerino, Italy 

(Dated: February 1, 2008) 

Abstract 

We demonstrate that spin chains are experimentally feasible using electrons confined in micro- 
Penning traps, supplemented with local magnetic field gradients. The resulting Heisenb erg-like 
system is characterized by coupling strengths showing a dipolar decay. These spin chains can be 
used as a channel for short distance quantum communication. Our scheme offers high accuracy in 
reproducing an effective spin chain with relatively large transmission rate. 

PACS numbers: 03.65.-w, 03.67.-a, 03.67.Hk, 03.67.Lx 
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I. INTRODUCTION 



Recently much theoretical research work has been focused on the possibility to use sys- 
tems of spins, coupled by ferromagnetic Heisenberg interactions and arranged along chain 
structures, for transferring quantum information. The remarkable property of these systems 
is the capability of transmitting the qubit state along the chain with fidelity exceeding the 
classical threshold and by means only of their free dynamical evolution. After the seminal 
paper by Bose in which the potentialities of the so-called spin chains have been shown, 
several strategies have been proposed to increase the transmission fidelity and even to 
achieve, under appropriate conditions, perfect state transfer jsl, Q, (|. All these proposals 
refer to ideal spin chains in which only nearest-neighbor couplings are present. However, also 
the more realistic case of longrange couplings, in particular magnetic dipole like couplings, 
has been studied 0, Qj- In [7] it has been shown that perfect state transfer or, at least, 
high transmission fidelity can be obtained by appropriately choosing the system parame- 
ters, such as local magnetic fields and inter-spin distances. Moreover, even when no site 
specific locally-tunable fields are allowed, spin chains with dipolar couplings often perform 
better, in terms of transmission fidelity, than their nearest-neighbor coupled counterpart . 
Hence, from these theoretical predictions, we expect that spin chains, also in the case of long 
range interactions, may represent a very promising system to achieve high fidelity quantum 
information transfer without requiring experimentally demanding gating operations. 

In this paper we demonstrate that a linear array of electrons, confined in micro-Penning 
traps, can implement an effective spin chain with magnetic dipole like spin coupling. The 
same system consisting of trapped electrons in vacuum has been already pro pos ed as a valid 
and competitive candidate for universal quantum information processing [9I. Iiol. 11]. In this 
respect, the possibility of reliably transmitting the qubit state between different quantum 
registers, without applying gate operations, is highly desirable. In fact, the use of a quantum 
channel to transfer a qubit state in a quantum processor can be a valuable alternative to 
the repeated application of swapping gates. 

We have already proved in [ll| that the addition of a magnetic field gradient, together 
with the Coulomb interaction between the particles, allows to obtain an effective nuclear 
magnetic resonance (NMR)-like coupling between the spins of the confined electrons. Here 
we generalize this approach to encompass a variety of trap set up, also in connection with 
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novel geometries of Penning traps 12|, [l3| . Indeed, by further investigating the interaction 



between the internal (spin) and external (motional) degrees of freedom of each particle, in- 
troduced by the applied local magnetic field gradient, we can mimic more general systems, 
with Heisenberg ferromagnetic or antiferromagnetic Hamiltonian. This fact opens up the 
possibility to simulate quantum spin systems with tunable interactions, thanks to the experi- 
mental control over the different trap parameters. The ultimate goal may be the observation 
of quantum phase transitions, as proposed with trapped ions controlled by laser beams Jjj]. 

In our proposal, we consider a linear array of electrons with inter-particle distance ranging 
from few microns to 50 fim. We provide an analytical expression for the spin-spin coupling 
strength, which shows a dipolar decay law. We estimate the value of the spin-spin coupling, 
for different ranges of the system characteristic frequencies as well as of the intensity of the 
magnetic gradient, with the aim of optimizing the transfer time of our quantum channel. 
Furthermore, we evaluate the fidelity of our system in reproducing an effective spin chain 
according to the Heisenberg model. In particular, we calculate the probability to obtain 
a perfect spin state transfer in a chain consisting of just two electrons. This probability, 
equal to one in an ideal spin chain l], in our system is less than one owing to the effects 
resulting from the interplay between the internal and the external degrees of freedom of the 
trapped particles. However, by an appropriate choice of the system parameters, especially 
the frequency and the amplitude of the spatial motions, we can obtain high fidelities and, 
at the same time, sufficiently large values of the spin-spin coupling. The electron trapping 
arrangement offers also the possibility to apply arbitrary site-specific changes in the system 
parameters in order to maximize, as outlined in 3, [8] , the efficiency of the quantum channel. 
Our theoretical predictions suggest that a linear array of electrons is suitable to implement 
a spin chain with the present technology. 

The manuscript is organized as follows. In Sec. [Til we describe the system and how the 
local magnetic field gradient couples the electron spin to the motional degrees of freedom. 
This coupling, mediated by the Coulomb interaction between charged particles, results in 
an effective Heisenberg-like Hamiltonian (Sec. IIIII) . In Sec. IIVI we estimate the fidelity and 
the efficiency of our system as a channel for quantum information transmission. Finally, the 
results of our analysis are summarized and discussed in Sec. |V] The more technical details, 
concerning the derivation of the fidelity, are presented in Appendix [A] 
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FIG. 1: Schematic drawing illustrating two different geometries for a linear array of micro- Penning 
traps. The traps are represented by sketching the electrostatic potential along the z axis, (a) The 
electrons are aligned along the z axis, parallel to the confining magnetic field; (b) the electrons are 
aligned along the x axis, orthogonal to the confining magnetic field. 

II. A LINEAR ARRAY OF TRAPPED ELECTRONS 



We consider a linear array of N electrons in micro- Penning traps [15]. According to the 
different geometry of the electrode arrangement, the micro-trap array can be either parallel 
to the direction of the confining magnetic field, i.e. along the z axis as shown in Fig. QJa), 
or orthogonal to this field, for example, along the shown in Fig. [T^b). To confine 

electrons in an array along the z direction we can use a closed cylindrical electrode structure 



9|, [10( or an open wire arrangement 



12] . This latter structure can also accommodate the 



electrons in an array aligned along the x axis. An orderly set of micro-traps, orthogonal to 
the trapping magnetic field, can be likewise realized by means of a planar electrode system 
[13I ] . As we will see, the different orientation of the linear array of particle affects the form of 
the resulting interaction Hamiltonian. Hence, we firstly derive the expression of the effective 
Hamiltonian in the case of micro-traps aligned along the z axis. Then we will show how this 
expression modifies in the case of an array directed along the x axis. 

The Hamiltonian of a system of N electrons confined in an array of Penning traps can 
be written as 

N 

H = T,Hr + T,H^ (1) 



i=l 



4 



where 



rNC (Pi - eA i) , deft 



H™ = V , + eVi - f^r< ■ B, (2) 

2m. 4m P 



represents the single electron dynamics inside each trap and 

Hfi = r ( 3 ) 

J 47re |ri — r j | 

describes the Coulomb interaction between electrons i and j. In Eqs. (J2J) and fl3j) m e , e, 
g, and cr^ are, respectively, the electron mass, charge, giromagnetic factor, and Pauli spin 
operators. As shown in Fig. Ufa), we assume that the micro-traps are aligned along the z 
axis and that z^o is the position of the center of the i-th trap. The electrostatic potential 

VtwH^V ^-vr-G + M 2 (4) 

is the usual quadrupole potential of a Penning trap, where V is the applied potential dif- 
ference between the trap electrodes and I is a characteristic trap length. The magnetic 
field 

Bi = -- (xi i + ya) + [.Bq + b(zi - z ifi )] k (5) 

is the sum of the uniform magnetic field B k, providing the radial confinement, with a local 
linear magnetic gradient around the i-th trap. The associated vector potential 

Ai = i [B + b(zi - z i)Q )] (-yii + Xij) (6) 



preserves the cylindrical symmetry of the unperturbed trapping field. 



llj , the spatial part of H, 
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Following an approach similar to the one described in Ref. 
can be recast in the form 

H- ext) ~ -Uuj m a} mi a m ^ + huj c a ] ci a c ^ + %LO z a\^a z>i 

+ fuj z e(a Z A + al A ) (—al ni a m ,i + al^cA , (7) 

' \U) C ' / 
where the annihilation operators a m ^ a Cj ;, a 2j ; refer, respectively, to the magnetron, 

cyclotron and axial oscillators of the i-th electron. The frequencies of these oscillators 

^ - 1^, (9) 
m e 

u - = (10) 



depend on the applied external fields and on the trap size. They build up a well defined 
hierarchy with uj m <C uj z <C u c . Indeed, we exploit this fact together with the assumption of 
a weak magnetic gradient, such that b\z{ — z^/Bq <C 1, to derive the Hamiltonian ([7j). The 
dimensionless parameter 

\e\b / h \e\bAz . . 

e = \ = , (11) 

m e ui z V 2m e u) z m e u z 

with Az being the ground state amplitude of the axial oscillator, represents the coupling, 
due to the magnetic gradient, between the axial motion and the radial degrees of freedom. 
In a similar way the magnetic gradient introduces also an interaction between the spatial 
and the spin motion. This coupling becomes evident by considering the spin part of the 
Hamiltonian Eq. (j2J) 

Tr(spin) 9 e ^ -r> 

gh , gh\e\b ,, . gh\e\b . _ v . 

= + 4m* ( * ~ Zifi) ~ 8^r {a ' Xi + GiVi) (12) 



which, in terms of the ladder operators, becomes 111 ] 



A 4 4 v LU r 



(13) 



where u c = Ju% — 2uj% is essentially a modified cyclotron frequency due to the insertion of 
the quadrupole potential. In deriving Eq. ( TTBl) we defined the operators = (of ±iof )/2 
and the spin precession frequency uj s = (g/2)u c . 

Hence, the Hamiltonian, Eq. (j2j), of a single electron can be written as 

h nc ^ -ft Wm aj fhi o mii + hu c a ] ci a c ,i + frc^a^a^ + ^ w »°f 

+ hu z e(a Z)i + a; J ^ H + 

V 4 



—ehu), 
4 



crJ +) (oc,i + aj, <) + o-J } (aj, + a m> j)l . (14) 



We now assume that, for each particle of the array, the cyclotron oscillator is in the ground 
state and the magnetron oscillator is in a thermal state with an average excitation number 



/ -C uj c /uj m . We recall that the ground state cooling of the cyclotron motion 
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or electrons 
and ions 
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181 ] have been experimentally obtained. Under these conditions, we can neglect in Eq. (JT 



the coupling between the axial oscillator and the radial motion. We can further simplify 



Eq. (j!4p by means of the rotating wave approximation (RWA). Indeed, terms like a\ a) mi 
are rotating at a frequency u s — u m much larger than the anomaly frequency u a = u s — u c 
typical of terms like <r\ a c> i and, therefore, are negligible in RWA. Hence, the Hamiltonian 
Eq. (Tl4"l) becomes 

h 

h nc ^ -7iu m a} mi a m ^ + Tiu c a ] ci a c ^ + hu z al^a z ^ + t^s^I 

+ -ehu z (a z4 + a\ ) of - -etiw g J-^- (o"| +) a C)i + crj _) a^) . (15) 
4 v 7 4 y u; c v ' 

We see that the applied magnetic field gradient couples the different electron spin compo- 
nents to the axial and to the cyclotron oscillators. 

Let us consider the part of the Hamiltonian Eq. ([3]) describing the Coulomb interaction 
between two electrons trapped at the sites i and j. If the oscillation amplitude of the two 
electrons is much smaller than the inter-trap distance dij = \zi t o — Zj t o\, we can expand the 
interaction Hamiltonian in a power series and retain the terms up to the second order 

flg - - V + Wi - %) 2 - 0* - - (* - Vif] , (16) 

where Zi = z—Zifi. The Coulomb interaction produces three effects on the electron dynamics: 
(i) a displacement of the equilibrium position along the z axis, (ii) a shift of the axial 
resonance frequency, and (iii) a coupling between the motional degrees of freedom of different 
particles. The first two effects are rather small and can be taken into account by redefining 
the trap center position and the corresponding axial frequency. Therefore, in the remainder 
of this section we focus on the coupled dynamics of the two electrons 
c e 2 „ 

where the coupling strength 

, = ^ \_* (±*) 2 (18) 

l ' J 8ne Q m e u} z di tj hAireodij \d itj J 

amounts to the Coulomb energy times the square of the ratio between the axial ground 
state amplitude and the inter-particle distance. We have observed that each oscillator, 



axial, cyclotron and magnetron, is characterized by a typical resonance frequency. As a 
consequence, the coupling introduced by the Coulomb interaction between the degrees of 
freedom of different electrons is effective only for almost resonant oscillators. Therefore, in 
Eq. (|T7j) the terms that couple the cyclotron and magnetron motion of the two particles give 
negligible effects. Furthermore, we are not interested in the coupling between the magnetron 
motion of different electrons, since this mode is essentially decoupled from the other degrees 
of freedom. Hence, disregarding the magnetron motion, the part of the system Hamiltonian 
describing the Coulomb repulsion between electrons i and j reduces to 



Hfj ~ -2h£ id (a Z j + al ti )(a z j + a\ d ) + 2h^~ (a c<i a\j + 4,i a <y) • (19) 



In the case of a linear array of electrons, trapped in a direction orthogonal to the magnetic 
field, i.e. along the shown in Fig. \JJh), we can derive a similar expression for the 

Coulomb coupling 



Hfj ^ Hi,j( a z,i + a l,i)( a z,j + aL) - \a>c,iaL + alifld 

UJ C 

+ 3(a Cti a Ctj + a^ajj)] . (20) 



We emphasize that in the case of Eq. ( 1191) . referring to the vertical array of traps, the coupling 
between the cyclotron oscillators of different electrons represent a swapping of excitations, 
which basically conserves energy. The only terms that survive involve the creation of a 
quantum of excitation at the site j at the expense of the destruction of a quantum of 
excitation at the site % and viceversa. In the case of an horizontal arrangement of traps, 
Eq. (|20l) . this is, in general, no longer true. Indeed, even though the leading terms are 
preserving the energy of the two coupled cyclotron oscillators, we also note the presence 
of rapidly rotating terms which involve the simultaneous creation and annihilation of two 
excitations. However if t; id (uj z /oj c ) <C uj c the effects of these rapidly rotating terms are 
negligible (RWA) and the Hamiltonian (120]) becomes 



H ?,j - HiA a z,i + + a l,j) ~ ( a c,i a l,j + a U a c,j) ■ (21) 

ui c 



We also note that Eqs. ffl9|) and ([2T|) exhibit alternating signs in front of the coupling 
terms. As we will see in the next section, this results in a different kind, ferromagnetic or 
antiferromagnetic, of the effective spin-spin interaction. 
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III. EFFECTIVE SPIN-SPIN COUPLING 



In the previous section we have seen that the magnetic gradient induces, for each particle 
of the array, a coupling between the spatial and the spin motions. This interaction, mediated 
by the Coulomb repulsion between the electrons, gives rise to an effective spin-spin coupling 
between different particles . This effect can be brought t o lig ht by making an appropriate 



unitary transformation on the Hamiltonian of the system |19j. We seek a transformation 
that formally removes, in the single particle Hamiltonian, the coupling between the internal 
and the external degrees of freedom of each electron. Hence, we transform the Hamiltonian, 
Eq. p, as H' = e s He~ s with 



S 
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N 

£ 

i=l 



+ 



UJ, UJ 



UJ n V UJ, 



^z ( 



°i a cA - °i a c, 



(22) 



where uo a = uo s — uj c is the anomaly frequency. The unitary transformation changes the 
operators, to the first order in e, in the following way 



g uj z lo z (_) 



-e 



4 UJ n 



UJ r 



°i + n £ — \— K a c tl 
2 wju c v 



U l "'C.l / 1 



(+) 



(J 



(+) , 9 c (+) ( \ 



9 u z uj z z f 



0" a' 



4 cu a y cu c 

To derive the expressions above we made use of the expansion 

e vA Be - V A = B + B] + t [A [A B]] + t [A [A [A B]]] + 



(23) 
(24) 

(25) 

(26) 



(27) 



where A and B are two noncommuting operators and rj is a parameter. 

The single electron part, Eq. ( fl5l) . of the system Hamiltonian can be written, after applying 

the unitary transformation, as 



H 



'NC 



~ -hu m aL M md + huj c a J ri a c ,i + huj z a\ i a z , i + '-uj s cr, 



h 



J m u 'in,i u "m,i 



z^z i Uj z,% 



(28) 



where we have neglected second and higher order terms in e, which in the cases relevant to 
the present analysis is of the order of 10~ 2 . Nevertheless, these extra terms are derived in 
Appendix [A] and their influence on the performances of the system is discussed in Sec. IIV1 
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Let us now turn to the Coulomb part of the system Hamiltonian. The first term in 
Eq. ffl9|) becomes 



( a z ,i + al A - 9 -ea^j (a z>j + a) z>j - 9 -ea] ) . (29) 



Expression (129]) contains a term which represents an effective spin-spin coup 



ing between 



different electrons in the array. This effect was already pointed out in Ref. [111 ]. Moreover, 
we note that the unitary transformation enforces a coupling between the axial motion of the 
j-th electron and the spin of the i-th electron. This effect is smaller of a factor £i,j/ui z <C 1 
than the corresponding coupling [see Eq. (|T5l) ] between internal (spin) and external (axial 
motion) degrees of freedom of the same particle. The error introduced by neglecting these 
terms is estimated in the Appendix. 

The remaining term in Hamiltonian ( fl9|) transforms into 



u z / g uj z \u) z (-)W + 9 u z u z (+) 
uj c \ ' 4 u a V uj c \ 4 u a V uj c , 



+ 2h& tj — a\i - -e — J —<r\ a c>j - -re— J —a) . 30 
oj c \ ' 4 uj a y uj c J \ 4 u a V / 

From Eq. (130|) we see that the unitary transformation produces the term 



which represents a direct coupling between the spin motion of different particles. Also in 
this case, there are additional terms in expression fl30l) . that couple the spin of an electron 
to the cyclotron motion of the other electrons in the chain. In comparison with the spin- 
cyclotron interaction for the same particle [see Eq. ( fl5i) ]. this coupling is reduced of a factor 
£ijU z /oj c u a , which, in the range of the parameters considered here, is typically much less 
than one. For an estimate of the errors introduced by these terms we refer to the Appendix. 

Hence, summarizing the results of our derivation, we have an effective spin-spin coupling 
between the electrons with the spatial dynamics substantially decoupled from the spin dy- 
namics. Consequently the spin part of the system Hamiltonian can be written, in the case 
of a linear array of electrons along the z axis, as 

//: - i \"s°t - ^E(2^H - ./;»;' - J*y^i (32) 

i=l i>j 
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where 

J z - ( 9 )\ c 2 ( 9 Y ^ (33) 

In Eq. (1531) we used the relations u a ~ lCT 3 u; c and cD c ~ uj c . We obtain a spin-spin interaction 
that is antiferromagnetic (ferromagnetic) if it is transmitted by the cyclotron (axial) motion. 

The situation is completely different when the linear array of electrons is aligned along 
the x axis 

N * * N 

H'.-Y, \"s°t • tX ( 2J Zi°W ~ J i»* ~ J^laf) . (35) 

i=l i<j 

In this case, the sign of the Heisenberg like coupling is reversed. The ferromagnetic (antifer- 
romagnetic) interaction is associated to the cyclotron (axial) motion. Similar results were 
also found in the case of ions, in a Paul trap, driven by six counterpropagating laser beams 

Q. 



IV. A CHANNEL FOR QUANTUM COMMUNICATION 

The Hamiltonians fl32|) and fl35l) describe a system of N spins coupled through Heisenberg- 
like interactions. These Hamiltonians can transmit an unknown spin state, from the electron 
placed at one end of the linear array, to the electron placed at the other end of the array. 
The remarkable fact is that this quantum information transfer is realized only by means of 
the free dynamical evolution of the system, without requiring any external action by the 
experimenter during the transfer. 

Therefore, let us analyze the potentialities of our system as a quantum communication 
channel. In our scheme, the dependence of the spin-spin coupling strength on the system 
parameters is shown in Eqs. ( J33l) and ( |34l) . In particular, «/?-, are proportional to l/d%j, 
that is they decrease with the distance between the particles % and j according to the dipolar 
decay law. They also depend on the applied magnetic field gradient and on the characteristic 
frequencies of the electron motion. More specifically, the value of J^J (J^-) depends on the 
cyclotron (axial) frequency uj c (uj z ). As a consequence of this fact we can neglect JfJ with 
respect to J? - when the value of the ratio uj c /u z is sufficiently large, as in the case considered 
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in [11] . Differently, in this paper, we choose smaller values for the ratio uo c /uj z (generally 
about 20 or less), so that J?V is of the same order of magnitude of J? - or even larger. Indeed, 
one can easily check, from Eqs. (I3"3"j) and that when uj c /uj z — 18.8 it is possible to obtain 
an isotropic Heisenberg-like interaction with 2J^- = JfJ. 

Generally the time required to transfer a qubit, encoded in the spin state, along a Heisen- 
berg chain depends on the values of J^J, so that the larger the value of J^J the faster the 
transfer. Indeed, the state transfer time t ex in a Heisenberg chain, consisting of just two 



spins, is equal to 



*« = 4J^- (36) 



We assume that the particles in our linear array are equally spaced with d = d^+i and J xy = 
Jii+i- From Eq. fl34l) we see that J x?/ oc b 2 /(u^d 3 ). Hence, to speed up the transfer process 
we have to miniaturize the system, to increase the strength of the magnetic field gradient 
and to reduce the cyclotron frequency. However, the value of the cyclotron frequency cu c , 
depending on the confining magnetic field, cannot be decreased at will, since it should be 
sufficiently large to cool the cyclotron motion to its ground state. For example, at the trap 
temperature of 80 mK Jj]] it is sufficient a cyclotron frequency of the order of 10 GHz. The 
inter-particle distance d depends on the level of miniaturization of the trap. We consider 
d varying from few microns to 50 /im. Finally stronger local magnetic gradients are, in 
general, achievable by reducing the micro-trap size. 

Essentially, the effective Heisenberg-like Hamiltonians, Eqs. ( 1321) and ( 1351) . have been 
obtained by taking two steps: we applied an appropriate unitary transformation and then 
disregarded the residual coupling between the different degrees of freedom. Both these steps, 
in general, introduce errors which reduce the accuracy of our system in reproducing an array 
of particles interacting according to the Heisenberg model. In particular, we neglected terms 
representing a residual coupling between the spin and the motional degrees of freedom, as 
well as between the different spatial oscillators. In Appendix [Aj we analyse in detail the 
role of each of these terms. Here, we only present the most relevant part of this interaction 



N 



H r ~ e 2 ^2 



Ji t / uj 2 z \ t 



(37) 



a a 9 m.i ""171,1 ' \ A x I "arc,! 

_4:UJ c UJ a AU Z C \4:UJ c UJ a J 

which affects the spin frequency, introducing a dependence on the cyclotron and magnetron 
motion. As a consequence each particle acquires a different spin frequency with a finite 
linewidth due to the thermal state of the motional degrees of freedom. 
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In order to know how precisely our model can simulate an ideal Heisenberg system we 
introduce the system fidelity 

F={i> f \Tr ext {p(t)W f ), (38) 

where 

|V/> = e-*"*U>, (39) 

with |^o) being the initial state of the spin chain and H s is the Heisenberg Hamiltonian 
Eq. (1321) . The operator p(t) in Eq. (1381) represents the density operator of the electron 
chain, including the motional degrees of freedom, evolved at the time t according to the 
full Hamiltonian of the system Eq. ([1]). We also assume that initially the axial, cyclotron, 
and magnetron motions are prepared in thermal mixtures with, respectively, an average 
excitation number k, n and I. The reduced density operator, describing the spin state, is 
then obtained by tracing over the spatial modes of the electrons. 

The system fidelity can be analytically calculated. The details are provided in Ap- 
pendix |A] In general, the fidelity can be written as 



1-S r - e 2 Ss 



(40) 



where £ r and £5 represent, respectively, the errors due to the residual coupling, Eq. fl37l) . 
and to the canonical transformation. In the simplest case of just two electrons, we find 

'5 s (ni,li,n2,hy 



£ r = 1 - E E PnMPMPnMPiih) 

rai,Zi n2,h 



4Jxy 



(41) 



with Pfn(m), Eq. ( 1A7I) . being the occupation probability for the mth Fock state, 

Fd(() = l 1 1 - + 



cos(f)sin(f v ^+C 7 ) , sin 2 (fv / T+C 7 ) 



and 



E 2 LJ 2 



n 2 - n lj 



2u% 



1 + C 2 



{k - k) 



(42) 



(43) 



being the detuning between the two spin frequencies. The fidelity decreases because of this 
finite detuning, which is determined by the thermal state of the cyclotron and magnetron 
oscillators. Indeed, the error, Eq. (141j) . vanishes in the ideal case of zero detuning 5 S = 0. 
Also the error due to the canonical transformation 



£ s = l -(2k+l) 



UJ': 



(5n + 1) + 



{51 + A) 



(44) 
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A 


B 


d fum) 


50 30 


10 3 


10 3 


lo z /2tt(MRz) 


490 490 


490 1200 


730 4500 


6(T/m) 


350 600 


1800 1800 


1100 1100 


I 


0.01 0.1 


2 50 


0.15 1 


J x y (kHz) 


0.01 0.14 


35 1300 


2.5 100 



TABLE I: Table showing the values of the axial frequency uj z /2tt, the magnetic gradient b, the 
average magnetron excitation number I and the coupling strength J xy for different choices of the 
nearest neighbor distance d. In case A (B) we have T = 0.99 (J- = 0.999) and uj c /2tt = 8 GHz 
(lo c /2tt = 11 GHz). We suppose that the axial and cyclotron motion are thermalized with the trap 
environment at the temperature of 80 mK. 

becomes larger when the electron motion is relatively hot. From Eq. (I44p . we see that this 
error is proportional to the average excitation numbers k, n, and I. Therefore, to increase 
the system fidelity it is essential to cool, possibly to the ground state, the electron motion. 
This comes automatically for the cyclotron oscillator, when the trap is at a temperature 
below 1 K, whereas the cooling of the axial and magnetron oscillators requires appropriate 
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18| 



techniques 

We present a number of cases in Tabled when the fidelity approaches the value one. We 
see that, for T = 0.99 {T = 0.999) and the inter-particle distance d ranging from 50 fim 
(10 /im) to few microns, we have a coupling constant J xy in the range 10 Hz 4- 1.3 MHz 
(2.5 kHz 100 kHz). For example, in the case of d = 10 jum we obtain J xy = 35 kHz by 
taking a cyclotron frequency uj c /2tt = 8 GHz, an axial frequency u> z /27r = 490 MHz, and a 
magnetic gradient b = 1800 T/m. 

We also recall that the decoherence time of the spin state as well as the heating time 
of the spatial motions, estimated according to the model described in jiol . I21I . 22], is much 
longer than the transfer time t ex . This remains true also for moderate values of the coupling 
strength J xy , thus allowing the transmission of the qubit state across the chain within the 
decoherence time of the system. 

Finally we note that our system offers the possibility, in principle, to apply arbitrary site 
specific changes to its parameters, such as the inter-particle distance, the magnetic gradient 
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strength, and the magnetic field magnitude. Hence, as suggested in [7|, |8j, by means of these 
local modifications one can optimize the transmission rate and the fidelity of the chain. 



V. CONCLUSIONS 

In this paper we presented a scheme for implementing a spin chain with long range in- 
teractions by means of a linear array of electrons confined in micro-Penning traps. Both 
antiferromagnetic and ferromagnetic Heisenberg-like systems can be realized using a local 
magnetic field gradient, mediated by the electrostatic interaction between the trapped par- 
ticles. In particular, we derived an analytical formula for the strength of the spin-spin 
coupling, which determines the transmission rate of the channel, as a function of the rel- 
evant system parameters like the inter-particle distance, the cyclotron frequency, and the 
value of the applied magnetic gradient. In our analysis we also estimated the fidelity of 
the system in reproducing an effective Heisenberg chain, by taking into account the effects 
produced by the coupling between the different degrees of freedom of the particles. We 
found that the fidelity depends on the frequency and the amplitude of the spatial motion of 
the particles. In general, higher values of the fidelity are obtained for smaller values of the 
spatial motion amplitudes and for larger values of the detuning between the characteristic 
trapping frequencies. The numerical estimates, calculated for an inter-particle distance d 
varying from 50 /mi to few microns, give a spin-spin coupling strength J xy in the range 10 Hz 
-=-1.3 MHz with a fidelity of 99%. Even in the case of a relatively weak coupling constant, 
the transmission of the qubit state from one end to the other of the chain takes place well 
within the decoherence time of the system. Moreover, the geometry of the system offers the 
possibility to apply arbitrary site-specific changes of its parameters in order to optimize the 
transmission rate and the fidelity of the quantum channel. 

In conclusion, an array of electrons confined in micro-Penning traps lends itself to im- 
plement, within the reach of current technology, quantum channels with high accuracy and 
sufficiently large transmission rates. Furthermore, the versatility of our scheme allows one 
to simulate also more general spin systems, in one and two dimensions, thus paving the way 
towards the observation of quantum phase transitions. 
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APPENDIX A: FIDELITY 



In this appendix we provide a brief description of the approach adopted to estimate the 
fidelity, as defined in Eq. (|38p . Our starting point is the complete single electron Hamiltonian 
Eq. (1141) . In order to remove from this Hamiltonian, to the first order in e, the coupling 
between the different particle motions, we should apply a unitary transformation which takes 
into account also the magnetron oscillator 



TV 

S = eJ2 

i=l 

g uj 



- < | | \ / t _ \ , 9^ \^x_ (Jr) n \ 



l^z ( 



_ (+) t 



4 0J a V uj, 



°* ( ( 



* (Jr-) 



(Al) 



This unitary transformation represents a generalization of the transformation Eq. ( l22|) . since 
it encompasses all the degrees of freedom of the particles. 
From the definition of the fidelity, Eq. (1551) . we can write [ijj] 



T = (^ / |Tr ex 4e- 5 e-^'*e s p(0)e- 5 e^' t e i ]|^ / ) 



S„iH't„Si 



(A2) 



where 



with 



H = H ext + H s + H r , 



N 

H ext = (-^mfll^m,. + huJcdl ^i + %UJ z a\ t 



(A3) 



(A4) 



i=l 



being the Hamiltonian describing the uncoupled external dynamics of the particles. The spin 
Hamiltonian H s is given in Eq. fl32l) . whereas H r includes the residual coupling between the 
spin and the spatial degrees of freedom 
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We assume that initially the cyclotron, axial, and magnetron oscillators are in a thermal 
mixture, each one represented by the usual density operator 

+00 

Pth = Prn(m)\m)(m\, (A6) 

m=0 



with [23 1 



Pm(m) = (-± = ) ( (A7) 
VI + m) VI + mj 

being the occupation probability of the mth Fock state of a harmonic oscillator with average 
excitation number m. The initial spin state of the chain is 

|^o) = (cos 1 1 |)i + e* sin ~ | T)^ | |) 2 . . . | [} N . (A8) 

The information is stored in the state of the first qubit and should be transmitted to the 
opposite end of the chain, to the Nth spin. Therefore, the density operator of the system 
at time t = is 

p(0) = p^ l ^®p^ ext \ (A9) 
The ideal final state of the spin chain is represented by the state vector 

|V/) = exp(~ff a t)|V>b), (A10) 

which is obtained from the initial spin state, Eq. (1A8j) . when the system is described by the 
Heisenberg Hamiltonian H s , Eq. fl32|) . 

Now to calculate the value of the fidelity, we make an expansion of Eq. ( 1A2I) in powers 
of S and consider terms up to the second order in e 

T ~ (ApA- 1 ) + ^(ApA^S 2 + ApS^- 1 + AS^A' 1 + S 2 ApA~ l ) 

- (ApSA^S - ASpA- l S + ASpSA- 1 + SApA^S - SApSA' 1 + SASpA^lAll) 

where we defined A = exp[—(i/h)H't], p = p(0) and (...) = {ipf\Ti ext [. . The first 

order terms in S have been omitted since their contribution, after tracing over the spatial 
degrees of freedom, is zero. 

In the absence of the residual couplings, contained in the Hamiltonian H r , the spin chain 
evolution is unaffected by the thermal state of the motional degrees of freedom, because 
[H ext ,H s ] = 0. This leads to (ApA^ 1 ) = 1. The corrections to the fidelity come both from 
the presence of H r and from the canonical transformation, represented by the remaining ten 
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terms of Eq. (lAllO . In order to separate the two effects, we first evaluate the impact of the 
unitary transformation when H' ~ H ext + H s . This greatly simplifies the procedure and 
allows to achieve an analytical expression for the fidelity 



T ~ (ApA^ 1 ) - e 2 £ s , 



(A12) 



where 
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E 

i=i 
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(A13) 



(o; 8 — u; m ) 2 

The expectation values (. . .)o and (•••)/ are calculated, respectively, over the initial and final 
state of the spin chain. At the swapping time, when the state of the first spin has moved to 
the other end of the chain, 

N 



N 
i=l 



$>n/ = -(iV-l)-cos0, 



N 
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JV 



EKO/I 



i=l 
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TV 
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sin 6 1 . 



i=l 



(A14) 

(A15) 

(A16) 

(A17) 
(A18) 



Moreover, after averaging over all the initial states in the Bloch sphere, i.e. evaluating 
(l/47r) J, Jq W Ss sm6d6d<j), we obtain 



£ s ~\(2k + l) + \^ 
6 bu c 



UJ 



-§ (2n + 1 + 3(A^ - l)n) + 



^,1 



(u; s — Ct> m )' 



(2Z + l + 3(iV-l)(Z + l) 

(A19) 

where fc, n, and / denote, respectively, the average axial, cyclotron and magnetron excitation 
number. The expression (1A19I) gives the error due to the unitary transformation. 

Let us consider now the effects of the Hamiltonian H r , contained in the term (ApA^ 1 ) 
of Eq. (1A12I) . The residual couplings produce mainly two effects: they induce transitions 
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between the motional states of the electron and make the electron spin frequency depend 
on the state of the external degrees of freedom. Both these effects, as we will see, reduce 
the system fidelity. 

The probability to observe transitions between the states of the electron motion can be 
easily estimated using a perturbative approach. Indeed, the probability for the transition 
l^m) ~~ > iV'n) is n °t larger than roughly A\(^ n \AH\^ m )\ 2 / {Tiuj nm ) 2 , where AH is the per- 
turbation, uj nm the transition frequency and \ipi) the i-th eigenstate of the unperturbed 
Hamiltonian. In our case, the Hamiltonian H r plays the role of AH and the terms, respon- 
sible for the transitions between electronic states, are in the last three lines of Eq. ( 1A5I) . For 
example, the terms proportional to a z ti a^a C} i induce transitions between the eigenstates 
\n, k, I, I) and \n — 1, k — 1, 1, f) of the single electron Hamiltonian 



with probability of the order of e [u^/u c (u z — uj a ) 2 ]kn. A similar perturbative approach 
allows us to estimate also the transition probability due to the other terms of Eq. (1A5I) . 
These probabilities, for the terms involving couplings between different motions of the same 
particle, are proportional to e A /(Auj) 2 where Auo denotes the detuning between the electron 
oscillation frequencies. Hence, the error, in this case, is always negligible because is a 
correction of the fourth order in e and, moreover, the characteristic frequencies of the electron 
motion are quite different from each other. Very small errors are also produced by the 
terms in Eq. flA5j) involving the dynamics of different particles. In this case the transition 
probabilities are of the order of e 2 {£ij/uj z ) 2 and £ 2 (£i ) j/uj z ) 2 (ui z /uj c ) 3 (ui z /uj a ) 4: . Indeed, these 
values, for our choices of the system parameters, are negligible. 

In addition to the state transitions, the residual couplings enforce a dependence of the 
spin frequency on the state of the particle motion. The correction Au s to the spin frequency 



depends on the cyclotron and magnetron excitations. Indeed, the constant shift proportional 
to uj 2 /2u) c u) a equally affects all the spins in the chain and, therefore, does not introduce any 
detuning between the spin frequencies. This term can be taken into account by redefining 
the spin precession frequency uo s . On the contrary, the last two terms of Eq. (1A21|) introduce 
a detuning between the spin frequencies along the chain, since the cyclotron and magnetron 



H 



—fiu m a* m a m + hu c a) c a c + huo z a\a z + -uJ s a z 



(A20) 




(A21) 
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oscillators are in a thermal mixture with fluctuating excitation numbers n and I. As a 
consequence each spin in the chain acquires a different frequency depending on the thermal 
state of the electron motion. This leads, as we will show, to a reduction of the system 
fidelity 

For the sake of simplicity, we restrict our analysis to the case of just two electrons in the 
chain. The corresponding Hamiltonian reads 



H, tl = -u x a{ + -u 2 o z 2 + 2hJ xy (a{ +] '<7< _) + cr^a^) - 2hJ z a z a z 2 , 



(A22) 



2 1 2 

with uji = uj s + Au(rii, k). The unitary evolution of the system gives at the swapping time 
t ex = n/4J*v 
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where ( = 5 s /(U xy ) with 

8 8 {n u h, n 2 , l 2 ) =uj 2 -uji = e 2 u z 
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2u : 
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(A25) 



Hence, by using the relations (1A23j) and (1A24j) . we obtain the system fidelity 



^(0 



1 + cos(^)sin(| v ^+C 7 ) + sin^fy^+C 7 ) 



(A26) 



VTTC* 1 + C 2 

when the cyclotron and magnetron oscillators are in the states \rii, with % = 1, 2. Conse- 
quently, in the case of a thermal mixture, the expression for the fidelity becomes 

'5 s (n 1 ,l 1 ,n 2 ,ky 



^^EE Pn{n,)Pi{k)Pn{n2)Pi{k) 

ni,h n,2,h 



4Jxy 



(A27) 



We use this formula together with Eq. (1A19I) to numerically evaluate the fidelity of our 
system. 
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